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SUMMARY 



Cornification and epidermal barrier defects are associated with a number of clinically diverse skin disorders. However, a suitable in vitro 
model for studying normal barrier function and barrier defects is still lacking. Here, we demonstrate the generation of human epidermal 
equivalents (HEEs) from human embryonic stem cells (hESCs) and induced pluripotent stem cells (iPSCs). HEEs are structurally similar 
to native epidermis, with a functional permeability barrier. We exposed a pure population of hESC/iPSC-derived keratinocytes, whose 
transcriptome corresponds to the gene signature of normal primary human keratinocytes (NHKs), to a sequential high-to-low humidity 
environment in an air/liquid interface culture. The resulting HEEs had all of the cellular strata of the human epidermis, with skin barrier 
properties similar to those of normal skin. Such HEEs generated from disease-specific iPSCs will be an invaluable tool not only for 
dissecting molecular mechanisms that lead to epidermal barrier defects but also for drug development and screening. 



INTRODUCTION 

The major function of the skin is to form a permeability 
barrier between an organism and its external environment. 
Integrity and cohesion of the "bricks and mortar" structure 
of the outermost layer of epidermis, the stratum corneum 
(SC), is essential for terrestrial life, the survival of which ul- 
timately depends on maintenance of permeability barrier 
homeostasis (Elias, 1983; Kalinin et al., 2002; Segre 2003). 
The SC is composed of corneocytes ("bricks") encased in 
a lipid-rich extracellular matrix (ECM; "mortar") (Nemes 
and Steinert, 1999). 

Corneocytes, which are composed of keratin microfi- 
brils surrounded by cornified envelopes (CEs), are the 
final products of a linear keratinocyte differentiation 
pathway from a mitotically active stratum basale (SB) 
through a transcriptionally active stratum spinosum (SS) 
and a stratum granulosum (SG) that transitions into anu- 
cleate SC cells. The CEs are formed from precursor pro- 
teins directly beneath the plasma membrane. Following 
membrane disintegration, Ca 2+ influx activates transgluta- 
minase, which irreversibly crosslinks CE proteins around 
filaggrin-associated keratin filaments. These CEs combine 



with secreted and processed lipids to form a functional 
epidermal barrier. 

The lipids are packaged into lamellar bodies in the SG 
cells. The lamellar body, with all its contents, is secreted 
in response to barrier perturbation. The secreted lipids of 
the SC are processed into lamellar membranes (Figures 1A 
and IB). Lamellar bodies also contain proteolytic enzymes 
and antimicrobial peptides, which when secreted along 
with lipid contribute to permeability barrier function 
(reviewed in Elias, 2012). 

Ichthyoses and atopic dermatitis are due to polymor- 
phisms or mutations in genes that control keratinocyte 
differentiation, cornification, or lipid metabolism. Lhese 
changes lead to defects in epidermal permeability barrier 
homeostasis that range from mild (can be easily mistaken 
for normal dry skin) to life threatening (Segre 2006; Smith 
et al., 2006; Richard, 2004). Although ichthyosis-like dis- 
eases are also present in animals and can be mimicked 
in transgenic mouse models, there are still no suitable 
in vitro models for these diseases. 

Human epidermal equivalents (HEEs) are in-vitro-gener- 
ated 3D models that are widely used in experimental set- 
tings. However, their utility for addressing the mechanisms 



CrossMark 



Stem Cell Reports 

Epidermal Permeability Barrier 3D In Vitro Model 




Figure 1. Epidermal Permeability Barrier within the SC of Normal Human Skin 

(A) Schematic drawing of normal human skin, depicting the epidermal permeability barrier within the SC. 

(B-F) TEM demonstrates the components of a competent permeability barrier in normal human skin. Lipids are packaged into lamellar 
bodies (LBs) in the SG cells. 

(B) LBs (arrows) are present in the uppermost SG cells. These LBs are not discrete droplets, but rather are part of a tubular network 
extending from the Golgi. 

(C) LBs, with aLl their contents, are secreted into extracellular space (asterisks) between the uppermost SG cells. De, desmosome. 

(D) Close-up of LB morphology. 

(E) Lipid from LBs is processed into lipid bilayers (asterisk) by SC enzymes that require acidity, thus forming lipid structures that are 
impermeable to passage of water and ions. 

(F) Ceramide-based lipids bind to the CE (asterisk), forming the cornified lipid envelope (CLE). 
Scale bars are noted in each panel. 



of various skin disorders, or for drug development and 
testing, has been limited by the fact that previously engi- 
neered HEEs do not form a fully developed epidermal 
barrier. In spite of advances in HEE engineering, such as 
integrating melanocytes (Nissan et al., 2011), macrophages 
(Linde et al., 2012), or dermal fibroblasts (Itoh et al., 2013) 
and pluristratified epidermis, issues involving the genera- 
tion of a functional permeability barrier in vitro remain 
unresolved. Further, in vitro studies have been limited 
by the fact that only a limited number of HEEs can be 
generated from one sample of epidermis, and the primary 
keratinocytes generated from this sample may contain pre- 
viously unidentified polymorphisms in genes that modify 
epidermal growth, differentiation, or barrier development. 
HEEs generated from immortalized keratinocytes develop 
even less well than those generated from primary keratino- 
cytes (Gotz et al., 2012). Therefore, in order to develop a 
HEE model that can be produced in an unlimited number 
of genetically identical units, we turned to human embry- 
onic stem cells (hESCs) and induced pluripotent stem 
cells (iPSCs), primary cells that are capable of infinite 
proliferation and whose genetic footprint can be fully 
characterized. 

We demonstrate here that by following an integrated, 
highly reproducible four-step protocol, we were able to 
differentiate hESCs/iPSCs into keratinocytes with gene- 
expression profiles similar to those of normal human 



keratinocytes (NHKs, i.e., primary keratinocytes isolated 
from biopsies of normal healthy skin). These hESC/iPSC- 
derived keratinocytes were then used to generate HEEs 
in an air/liquid interface culture exposed to sequential a 
high-to-low humidity environment. HEEs generated from 
hESC/iPSC-derived keratinocytes have a functional perme- 
ability barrier similar to native human skin and are indis- 
tinguishable from HEEs generated from NHKs under the 
same condition. We are not aware of any other protocols 
resulting in this functional endpoint. 

RESULTS 

Reprogramming 

Since the first iPSCs were derived, different reprogramming 
approaches aimed at increasing their efficiency and mini- 
mizing the potential for DNA damage have been developed 
(reviewed in Bellin et al., 2012; Robinton and Daley 2012). 
We used a nonintegrating, nonviral, synthetic, modified 
mRNA system (Yakubov et al., 2010; Warren et al., 2010) 
to reprogram human newborn fibroblasts (Figure SI avail- 
able online). In our hands, the reprogramming efficiency 
was ~0.76% (Figures 2A and S2), which is lower than that 
reported by Warren et al. (2010). We chose two iPSC clones, 
iKCL004 and iKCLOll, that had a population doubling 
time (PDT) similar to that of KCL034, a clinical-grade 
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hESC line that we used as a control throughout all experi- 
ments. Both iKCL004 and iKCLOll had a methylation 
signature of PSCs, expressed pluripotency markers, and 
were able to differentiate into cells of all three germ layers 
in vivo (Figures 2B-2G and S3). Array comparative genomic 
hybridization (CGH) did not detect any copy number 
changes, using Promega male as a standard, over an 
extended period of cell culture (>3 months) under feeder- 
free conditions. 

Efficient and Reproducible Differentiation of hESCs 
and iPSCs into Keratinocytes 

To recapitulate epidermogenesis, in the first step of kerati- 
nocyte differentiation, which we termed "induction," we 
transferred a clinical-grade hESC line (KCL034) and two 
iPSC lines (1KCL004 and iKCLOll) from 5% O z to atmo- 
spheric 0 2 , and exposed them to 25 ng/ml BMP4 and 
1 u.M of the carboxylic acid form of vitamin A, all-fra«s 
retinoic acid (ATRA), for 7 days (Figure 3A). The lines did 
not respond similarly to induction. Both iPSC lines ex- 
pressed higher mRNA levels of the epidermal markers 
integrin (34 (ITGB4) and keratin 14 (KRT14), suggesting 
differentiation toward an epidermal lineage, whereas 
hESC line KCL034 seemed to remain "undecided," having 
similar mRNA levels of these two epidermal markers and 
a neural marker, oligodendrocyte lineage transcription 
factor 2 (0LIG2; Figure 3B). In the next step, termed "selec- 
tion," we used native 3D decellularized human dermal 
fibroblast (HDF) ECM as a growth-supporting platform. 
This was done to mimic basement membrane (BM), which 
is formed soon after epidermal induction during embr- 
yogenesis (Figure S4). A week later, we purified epidermal 
progenitor stem cells based on their preferential adherence 
to collagen IV (Murray et al., 1979; Bickenbach and Roop, 
1999). During the next "enrichment" step of our differenti- 
ation protocol, the purified cells gave rise to a homoge- 
neous population of keratin 14 (K14) + p63 + colonies, which 



was further amplified in the final "expansion" step of the 
protocol (Figure 3C). p63 is a transcription factor that plays 
a major role in maintaining the proliferative potential of 
epithelial progenitors (reviewed in Blanpain and Fuchs, 
2007). Whereas plating of cells on 3D HDF ECM rapidly 
increased p63 mRNA expression, pluripotent markers, 
such as the transcription factors Oct4 and nanog, steadily 
decreased during differentiation (Figure 3D). K18, which 
is expressed earliest during embryonic development, and 
with its binding partner K8 can serve as a marker of simple 
epithelium, was expressed at higher levels only during the 
induction and selection steps, suggesting that the cells sub- 
jected to our differentiation protocol were following a fate 
similar to that which occurs for epidermal morphogenesis 
during embryonic development. This four-step protocol 
was highly reproducible. KRT14 expression increased 
steadily during differentiation and there were no signifi- 
cant differences among the three lines in multiple rounds 
of the differentiation process (Figure 3E). The PD of 
hESC/IPSC-derived keratinocytes did not differ much 
from that of NHKs over a period of 5 weeks in culture 
(Figure 3F). 

hESC/iPSC-Derived Keratinocytes Are Very Similar to 
Primary NHKs 

Next, we compared enriched and expanded populations 
(T3) of hESC/iPSC-derived keratinocytes with NHKs from 
different donors in several different ways, and did not 
find any major discrepancies. Differences in the expression 
levels of KRT14 and p63 in multiple T3 populations of 
hESC/iPSC-derived keratinocytes (n = 9; three biological re- 
peats for each of three lines) and NHKs from three donors 
were not statistically significant (Figure 4A). Flow-cytome- 
try analysis of NHK and T3 populations of hESC/iPSC- 
derived keratinocytes using K14 and integrin f54 as markers 
did not identify any marked difference in the percentage 
of double-positive (DP) cells; 95% of NHK, 87.8% of 



Figure 2. Generation and Characterization of Normal Human iPSCs 

(A) Trypan blue staining of colonies in one well of a six-well dish. The average reprogramming efficiency from two independent experi- 
ments was ~0.76%. 

(B-D) DNA methylation analyses. Hierarchical clustering (B), methylation portrait (C), and principal component analysis (D) suggest that 
both iKCL004 and iKCLOll have the methylation signature of PSCs. BJ, parental HFFs. n = 3 biological replicates (rounds of differentiation) 
for each cell line. 

(D) Undifferentiated iKCL004 and iKCLOll cultured on HFF feeders and under feeder-free conditions display the typical morphology 
of PSCs. 

(E) Expression of pluripotency markers is confirmed by alkaline phosphatase (AP) activity assay and immunostaining for NANOG, 0CT4, 
TRA-160, and TRA-1-81 in iKCL004 and iKCLOll lines. 

(F) Differentiation of iKCL004 and iKCLOll into alL three germ layers in vivo is confirmed by detection of specific markers. Teratomas were 
encapsulated and did not invade surrounding tissues. ALL sections were stained with hematoxylin and eosin (H&E). Specific stains are 
brown (immunohistochemistry) or light blue (ALcian). Positive staining for mitochondrial complex IV type II confirms the human origin of 
the teratoma tissue. Germ layer markers: Alcian blue/periodic acid Schiff (PAS)-stained cartilage and desmin for mesoderm, filll tubulin 
and glial fibrillary acidic protein (GFAP) for ectoderm, and GATA4 and a-fetoprotein for endoderm. 

See also Figures S1-S3. 
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Figure 3. Efficient and Reproducible Dif- 
ferentiation of hESCs and iPSCs into Ker- 
atinocytes 

(A) Schematic of the differentiation proto- 
col. TO (day 0), Tl (day 7), T2 (day 14), and 
T3 (day 28) represent time points at which 
the cells were switched to a new culture 
condition and/or collected for analyses. 

(B) qPCR analysis measuring the expres- 
sion of differentiation markers (epidermal: 
ITGB4 and KRT14; neural: 0LIG2) suggests 
a difference in commitment toward epi- 
dermal lineage at time point Tl. n = 9 for 
each cell line; each of three rounds of 
differentiations (biological replicates) had 
three technical replicates. 

(C) Diversity is lost at later stages, and cell 
colonies at the enrichment step are uni- 
formly positive for keratinocyte markers 
K14 (red) and p63 (green). 

(D) qPCR analyses measuring expression of 
differentiation markers at TO, Tl, T2, and T3 
of the differentiation protocol. Expression 
of pluripotency markers 0CT4 and NANOG 
is diminished by the end of differentiation, 
whereas KRT18 and p63 mRNA levels reflect 
stages during epidermogenesis. n = 9 for 
each cell line; each of three rounds of 
differentiations (biological replicates) had 
three technical replicates. 

(E) Levels of KRT14 mRNA expressed during 
three rounds of differentiation at four 
different time points for each of lines 
demonstrate the reproducibility of the 
protocol. Absolute quantification was per- 
formed using the ddCt method. GOI ex- 
pression was normalized to the geometrical 
mean of two HK genes (GAPDH and EIF4A2). 
The SEM was calculated for mean expression 
in KCL034, iKCL004, and iKCLOII. n = 9 for 
each cell line; each of three biological 
replicates/rounds of differentiations had 
three technical replicates. 

(F) The proliferation rate of hESC/iPSC- 
derived keratinocytes shows no marked 

difference compared with freshly isolated primary NHKs over a period of 5 weeks. Statistical significance was calculated for each passage 
using two-way ANOVA with Sidak's multipLe comparison test. Data are presented as a blank adjusted mean with SD (n = 3 biologicaL 
replicates for each cell line). 
See also Figure S4. 
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KCL034 T3, 95.4% of iKCL004 T3, and 86% of iKCLOII T3 
were K14 + /integrin |34 + (Figures 4B and S5). Finally, we used 
the HumanHT-12 Expression BeadChip whole-genome 
expression array, which provides coverage for more than 
47,000 transcripts and known splice variants, to compare 
the transcriptomes of undifferentiated hESCs/iPSCs (TO), 
hESC/iPSC-derived keratinocytes (T3), and NHKs from 



two different donors. The TO and T3 transcriptomes were 
analyzed as biological replicates from three independent 
differentiation runs for each of three cell lines (KCL034, 
iKCL004, and iKCLOII). Principal component analysis 
identified an expected pattern: hESC/iPSC-derived kera- 
tinocytes at T3 were grouped together with NHKs and 
were dissimilar to undifferentiated hESCs/iPSCs at TO 
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Figure 4. Profiling of hESC/iPSC-Derived Keratinocytes 

(A) qPCR analyses measuring expression of KRT1 4 and p63 in hESC/iPSC-derived keratinocytes atT3 (n = 9; three rounds of differentiation 
for each of three lines) and NHKs (n = 2) found no significant difference between the two groups. Multiple t tests. Statistical significance 
was determined using the Holm-Sidak method. For KRT14, p = 0.120881; for p63, p = 0.155472. 

(B) Populations of hESC/iPSC-derived keratinocytes at T3 have a similar percentage of K14 + /integrin p4 + cells as NHKs. The percentage is 
an average from three rounds of differentiation for each of the three lines. At TO, no K14 + /integrin (34 + cells were detected. Fluorescence- 
activated cell sorting (FACS) images represent one of the three rounds. 

(C and D) Scatterplots showing the relative expression of all genes on the chip. 

(C) Global gene-expression profiles of NHKs from two donors (NHK1 and NHK2) shown as a scatterplot against each other. 

(legend continued on next page) 
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(Figure S6). As an average of three biological replicates, the 
transcriptome of iKCLO 11 -derived keratinocytes at T3 was 
closer to the transcriptome of NHKs than to keratinocytes 
derived from iKCL004 or KCL034. A scatter analysis of 
gene expression in NHKs from two different donors 
(NHK1 and NHK2) found that 939 transcripts (2.00%) 
differed by >2-fold (Figure 4C). An average of these two 
was then compared with an average of three biological 
replicates from each of the three lines at T3. The keratino- 
cytes derived from iKCLO 11 differed by >2-fold in only 
155 transcripts (0.33%) from NHKs by >2-fold, whereas 
iKCL004 (636 transcripts, 1.35%) and KCL034 (1,604 
transcripts, 3.41%) were more distinct (Figure 4D). Among 
all transcripts that were different from NHKs, only 122 
(0.26%) were common to all hESC/iPSC-derived keratino- 
cytes at T3 (Figure 4E). We profiled these 122 transcripts 
using the GeneGo database (Thomson Reuters) for 
Pathway Maps, Process Networks, Gene Ontology Net- 
works, and Disease Biomarkers list, and found that 34 
(0.07%) of them might have relevance to keratinocyte 
biology (Tables S1-S5, S6, and S7). 

We also looked at how reprogramming and differentia- 
tion changed human leukocyte antigen (HLA) gene expres- 
sion (Figure 4F). Although parental HDFs do express HLA 
class I, the reprogrammed clonal lines iKCL004 and 
iKCLOll do not, similarly to hESC line KCL034. Differenti- 
ation into keratinocytes changed the global pattern. 
Notably, HLA class I genes were turned on, implying that 
the recipients of epidermis generated from the hESC- 
derived keratinocytes as an off-shelf product would be in 
need of immunosuppressive therapy, whereas that may 
not be the case with iPSC-based autografts. 

hESC-derived keratinocytes exhibit an epidermal tran- 
scription program, and p63 is considered to be the main 
transcription regulator of epidermogenesis and essential 
for epidermal commitment of hESCs (Blanpain and Fuchs, 
2007; Aberdam et al., 2007; Medawar et al., 2008; Pozzi 
et al., 2009; Crum and McKeon 2010; Vanbokhoven 
et al., 2011). Therefore, we compared the expression 
profiles of 96 genes that, according to the NCI Pathway 



Interaction Database, constitute the p63 transcriptional 
network. Of these genes, only three (ADA [OMIM: 608958], 
SMARCD3 [OMIM: 601737], and YWHAQ [OMIM: 609009] 
had > 2-fold higher expression in NHK from two donors 
than in all biological replicates of hESC/iPSC-derived 
keratinocytes at T3. Similarly, only three of those 96 genes 
(BBC3 ox PUMA [OMIM: 605854], CDKN1A orp21 [OMIM: 
116899], and SP1 [OMIM: 189906]) had >2-fold higher 
expression in hESC/iPSC-derived keratinocytes at T3 than 
in NHKs from two donors. None of these six genes are 
known to be essential for epidermogenesis, although they 
may be indirectly involved. 

HEEs with a Functional Epidermal Permeability 
Barrier Generated In Vitro 

Molecular characterization suggested that our highly 
reproducible protocol for differentiating hESCs/iPSCs into 
keratinocytes yields a fairly pure population of cells that 
are very similar to NHK. To test the capability of hESC/ 
iPSC-derived keratinocytes to generate multilayered 
epithelium in vitro, we used the HEE model developed by 
R.S., A.C., D.C., M.H., L.C., S.D. Peckenpacker, K. Park, 
Y. Uchida, K.R. Feingold, P.M. Elias, D.I., and T.M.M. 
(unpublished data). The HEEs generated using this proto- 
col developed in vitro all of the cellular strata seen in the 
human epidermis. Since they also developed normal 
barrier properties in vitro, there was no need to implant 
the organotypic cultures into immunodeficient mice, 
which is a commonly used technique (Del Rio et al., 
2002). We assessed the integrity of the stratifying cultures 
by measuring the transepithelial electric resistance (TEER), 
which indicates changes in the transcellular and para- 
cellular permeability of epithelial/epidermal in vitro cell 
cultures. Within 2 weeks, the TEER in all cultures reached 
> 1,000 Q, indicating the formation of a permeability 
barrier and multilayered epithelia (Figures 5A and 5B). 

Differentiation markers were normally expressed and 
located (Figures 6A-6D). Since our system does not contain 
mixed fibroblast-keratinocytes cultures, the basement 
membrane zone (BMZ) generated by this system cannot 



(D) GLobaL gene-expression profiles of hESC/iPSC-derived keratinocytes at T3 (average of n = 3 bioLogical replicates/rounds of differen- 
tiation for each Line) shown in scatterpLots against NHKs (average of n = 2). Normalized Linear expression LeveLs for each probe are pLotted 
as XY coordinates for each sampLe. YeLLow Lines indicate boundaries of a 2-fold relative expression difference. 

(E) Venn diagram of the transcripts with a 2-fold relative expression difference from NHKs in KCL034-, iKCL004-, and iKCLOll-derived 
keratinocytes at T3 shows that 122 out of >47,000 were common to all three lines. 

(F) Heatmaps showing relative expression patterns across multiple sampLes for specific gene Lists. HLA genes were seLected by text 
searching for "HLA" within the "Gene SymboL" annotation, and p63-reLated genes were identified from the p63 transcriptionaL network 
in the NCI Pathway Interaction database. The coLors of the heatmaps indicate high (yeLLow) or Low (bLue) expression reLative to the 
average for that gene across the sampLes. The vaLues are the average of three independent runs for KCL034, iKCL004, and iKCLOll at TO 
and T3. Transcripts for NHKs from two different donors are shown independentLy on the HLA heatmap, and as an average on the p63 
transcriptionaL network heatmap. 

See aLso Figures S5 and S6, and TabLes S1-S5, S6, and S7. 
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Figure 5. Pluristratified HEEs Generated In Vitro from hESCs and iPSCs 

(A) The TEER during HEE formation over a period of 14 days reflects permeability barrier formation (n = 4 for NHK, n = 12 for each of lines; 
each point represents an average of measurements from three different spots). Error bars represent SD. Only 14-day-oLd HEE cultures that 
displayed a TEER of >1,200 ohm cm 2 were used for subsequent analyses unless indicated otherwise. 

(B) H&E staining demonstrates normal morphology with the presence of all epidermal layers regardless of the source of keratinocyte 
population. SB, stratum basaLe; SS, stratum spinosum. 



be equivalent to the BMZ of human skin. Nevertheless, in 
our 3D cultures, basal-layer keratinocytes were secreting 
ECM proteins, such as laminin, and assembling a BMZ- 
like structure (Figures 6E and 6F). 

Next, we examined the epidermal permeability barrier. 
Lipid secretion appeared normal. Lamellar bodies and 
extensive intercellular bilayers were detected in all cultures 
(Figure 7 A). The endoplasmatic reticulum (ER) Ca 2+ store, 
which is essential for normal keratinocyte signaling and 
differentiation (Callewaert et al., 2003; Celli et al., 2012), 
was detectable at day 8 of the culture (Figure 7B). Typically, 
the Ca 2+ gradient for normal epidermis has low levels in 
the SB and lowest levels in the SS, followed by an increase 
in extra- and intracellular Ca 2+ that peaks within the SG 
(Elias et al., 2002). In terminally differentiated air-lifted 
cultures at day 14, HEEs generated from NHKs and hESC 
(KCL034)-derived keratinocytes had the greatest density 
of Ca 2+ precipitates in the SG, whereas precipitates were 
somewhat more diffuse in HEEs generated from iPSC- 
derived keratinocytes (Figure 7C). However, the lanthanum 
perfusion assay, which is used to depict pathways of water 
movement through epidermis, demonstrated a functional 
permeability barrier in all cultures: the passage of elec- 
tron-dense lanthanum tracer was blocked by the epidermal 
lipid barrier at the SG/SC interface and was not detected 
in the SC (Figure 7D). 

DISCUSSION 

We have described a stepwise protocol for manufacturing 
in vitro HEEs, with all of the cellular strata and apparently 



normal barrier properties seen in the human epidermis, 
from hESCs and iPSCs. 

When manipulating the cellular fate of pluripotent cells 
in vitro, it is crucial to have a thorough understanding of 
the embryonic morphogenesis of the target tissue in order 
to achieve efficient production of the desired cell type (Hay 
et al., 2008; Cheng et al., 2012). Epidermal morphogenesis 
and homeostasis are relatively well described (Blanpain and 
Fuchs, 2009; Simpson et al., 2011). After gastrulation, the 
embryo surface emerges as a single layer of neuroectoderm, 
which will ultimately specify the nervous system and skin 
epithelium. This decision is guided by Wnt signaling, 
which blocks the ability of the ectoderm to respond to 
fibroblast growth factors (FGFs). In the absence of FGF 
signaling, the cells express bone morphogenetic proteins 
(BMPs) and become fated to develop into epidermis. In 
contrast, neural fate develops when, in the absence of a 
Wnt signal, the ectoderm is able to receive and translate 
activating cues by FGFs, which then attenuate BMP 
signaling through inhibitory cues (Chang and Hemmati- 
Brivanlou, 1998; Stern, 2005). Similarly to our protocol, 
most of the published protocols for differentiation of 
PSCs (either mouse or human) into keratinocytes follow 
this strategy (referenced in Figure S6). We introduced a 
step in which we purified epidermal progenitor stem cells 
based on their preferential adherence to collagen IV (Mur- 
ray et al., 1979; Bickenbach and Roop, 1999). In our hands, 
this step, which is a kind of physiological purification 
step that puts a minimal stress on the cells, resulted not 
only in a higher yield of hESC/iPSC-derived keratinocytes 
but also in a purified population that was highly homo- 
geneous and similar to NHKs (Figure 2). Homogeneous 
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populations of keratinocytes differentiated from disease- 
specific iPSCs would be advantageous for addressing the 
effects of such mutations on keratinocyte biology (Mildner 
et al., 2010; Simpson et al., 2010). 

Although a link between the level of ambient humidity 
and barrier development has been described in both ani- 
mals (Hanley et al., 1997; Komiives et al., 1999) and human 
premature infants (Agren et al., 2006), and then recapitu- 
lated in vitro with NHKs (Thakoersing et al., 2010), the 
exact molecular mechanism is still unclear. A limited num- 
ber of studies have addressed the relationship between 
environmental humidity and barrier function in postnatal 
epidermis (Takei et al., 2013; Katagiri et al., 2003; Sato et al., 
2002; Denda et al., 1998), and the results warrant further 
work. Regardless, the model that we describe here would 
have an unprecedented advantage for elucidating the 
molecular mechanisms of barrier development, perturba- 
tion, and recovery, as well as how mutations of genes 
involved in cornification and lipid metabolism affect 
permeability barrier homeostasis (Mildner et al., 2010; 
Simpson et al., 2010). 

Previous HEE models have developed some barrier prop- 
erties (Gschwandtner et al., 2013; Thakoersing et al., 2012; 
Ponec et al., 2001). The technique presented here is an 
advance from previously published HEEs because it is 
able to produce a HEE that not only shows morphologic 
changes associated with barrier establishment but also 
demonstrates the essential functional properties seen in 
skin with a functional barrier. The HEEs described above 
provide an effective block against water permeation, as 
evidenced by their ability to block lanthanum and form a 
calcium gradient. Further, this 3D model of HEEs with a 
functional permeability barrier is highly reproducible 
and has the potential to be easily scaled up and adapted 
to cGMP requirements for the use in regenerative and 
aesthetic medicine, as well as drug development. The 
recent (March 2013) European Union ban on the import 
and sale of cosmetics containing ingredients tested on 
animals will increase the demand for efficient and cost- 
effective alternative in vitro models that would not com- 
promise consumer safety. HEEs with a functional perme- 
ability barrier generated completely in vitro from either 
hESCs or iPSCs, as described here, might present the best 
available answer (Kuchler et al., 2013). 

EXPERIMENTAL PROCEDURES 

mRNA Synthesis and Modifications 

The intronless coding sequences (CDSs) of OCT4, SOX2, KLF4, 
c-MYC, and LIN28 were retrieved from the National Center for 
Biotechnology Information (NCBI). As previously described by 
Warren et al. (2010), a KOZAK translational initiation signal was 
added in the 5' UTR, whereas the 3' UTR was terminated with an 



alpha-globin-driven oligo(dT) sequence for templated addition 
of a polyA tail. The CDSs were then codon optimized with the 
use of GeneOptimizer expert software (GENEART/Life Technolo- 
gies). Codon optimization allows more efficient transcription 
by avoiding RNA secondary structures, adjusting GC content, 
and removing splice-sites. The gene encoding the destabilized 
enhanced GFP (d2EGFP) codon from pRVGP-d2EGFP vector 
(Addgene) and optimized OCT4, SOX2, KLF4, c-MYC, and LIN28 
CDSs were synthesized and subcloned (GENEART/Life Technolo- 
gies) into the pTNT plasmid vector, which was specifically de- 
signed for in vitro expression (Promega). Plasmids were linearized 
with BamHI. PCR templates were generated from the linearized 
plasmids using the same upstream primers for all six constructs 
and six different gene-specific, downstream T120-heeled reverse 
primers for addition of the poly(A) tail. Quality control of the 
PCR reaction was performed by capillary electrophoresis (Agilent 
Bioanalyzer). The templates were transcribed with T7 polymerase 
and further modified according to Warren et al. (2010). Modifica- 
tions included the complete substitution of 5-methylcytidine 
(5mC) for cytidine and the substitution of pseudouridine (psi) 
for uridine to increase mRNA viability and ectopic protein expres- 
sion. Furthermore, to promote efficient translation and increase 
mRNA half-life, a synthetic 5' guanine cap 3'-0-Me-m7G(5') 
ppp(5')G was added and 5'-triphosphates were eliminated by 
Antarctic phosphatase treatment. The DNA template was removed 
with DNase, and single-strand RNA (ssRNA) was purified with 
spin columns. Quality control of the purified ssRNA product was 
performed by capillary electrophoresis (Agilent Bioanalyzer; Fig- 
ure SI). An "mRNA cocktail" containing OCT4, SOX2, KLF4, 
C-MYC, LIN28, and c\2eGFP in molar ratios of 3 : 1 : 1 : 1 : 1 : 1 at a final 
concentration of 100 ng/ul was prepared (AmsBio), aliquoted for 
single use, and frozen in single-use aliquots. 

Reprogramming 

Mitotically inactivated foreskin fibroblasts (NuFF from donor 11; 
AmsBio) were used as feeder cells and plated at a density of 
2.5 x 10 5 /well of a six-well dish coated with CELLstart (Life 
Technologies). Human neonatal foreskin fibroblasts BJ (ATCC, 
CRL-2522) were plated onto feeders 24 hr later at a density of 
7,000/well. For each four wells to be transfected, 25 ul of RNAiMAX 
(Sigma), a lipofectamine transfection reagent, was diluted in 225 \A 
of Opti-MEM (Life Technologies) and mixed with 50 of mRNA 
cocktail diluted in 200 Opti-MEM. The complex was incubated 
at room temperature for 15 min before 120 ml was added dropwise 
into each of four wells to be transfected. The culture medium was 
refreshed 4 hr later with complete Pluriton medium (Stemgent) 
supplemented with 300 ng/ml B18R (eBioscience). Transfection 
was repeated for 17 consecutive days. On days 6-17, Pluriton was 
preconditioned with NuFF feeders for 24 hr and supplemented 
with 300 ng/ml B18R and 0.5 mM Na-butyrate. On days 18-21, 
potential iPSC colonies were left to expand in complete Pluriton 
with no other supplements. Ten colonies were picked based on 
morphology on day 21, transferred onto human foreskin fibroblast 
(HFF) feeders in KOSR-XF medium (Life Technologies), and 
expanded. Two clones with a similar growth rate, iKCL004 and 
iKCLOll, were further characterized. The cells were cultured at 
5% 0 2 , 37°C, 5% C0 2 throughout the reprogramming. 
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Characterization and Culture of PSCs 

hESC line KCL034 was derived and characterized as previously 
described (Jacquet et al., 2013; Stephenson et al., 2012). Expression 
of pluripotency markers and capability to differentiate into cells 
of three germ layers were confirmed for both iKCL004 and 
iKCLOll by methods that are routinely used for hESC lines in 
our derivation center. 

PSC lines were cultured under standard feeder-free conditions 
on BD Matrigel hESC-qualified matrix (Becton Dickinson) in 
mTESRl medium (StemCell Technologies) under hypoxic (5% 0 2 ) 
conditions. 

Methylation Assay 

Genomic DNA was isolated from BJ parental fibroblasts (negative 
control), KCL034 (positive control), iKCL004, and iKCLOll. Each 
sample of 500 ng DNA was bisulfite treated (EZ Zymo DNA Kit) 
according to the manufacturer's recommendations optimized 
for 450k array. Epigenome-wide methylation (>480,000 CpG sites) 
was analyzed using the Infinium Human Methylation 450 Bead- 
Chip Kit (Illumina) in biological triplicates for each line. GenoSplice 
Technology performed quality control, processing, and further 
analyses of the data. 

3D HDF ECM 

Preparation of decellularized 3D ECM was performed as previously 
described (Ilic et al., 2012). 

Differentiation Protocol 

To induce differentiation, undifferentiated PSCs were transferred 
into a 20% 0 2 atmosphere environment and treated with mTESRl 
basal media supplemented with 1 \iM ATRA (Sigma-Aldrich) and 
25 ng/ml BMP4 (R&D) for 7 days (Induction). To select for cells 
with early acquisition of ectodermal fate, the cells were harvested 
and replated onto freshly prepared 3D HDF ECM at a density of 
5-10 x 10 3 cells per cm 2 and grown in Dulbecco's modified Eagle's 
medium/Ham F12 (3:1; Life Technologies) supplemented with 
1 |xM ATRA and 25 ng/ml BMP4 for a further 7 days (Selection). 
To enrich for putative epidermal progenitors, rapid adhesion to 
type IV collagen-coated dishes was used, and the rapidly adhering 
cells were cultured in defined keratinocyte-SFM supplemented 
with 1 \iM ATRA for 7 days (Enrichment). After that, the cells were 
cultured in EpiLife medium (Life Technologies) for a further 
7 days (Expansion) before final harvest and analysis. 

Proliferation Assay 

All three lines and NHKs were seeded in six-well plates in triplicates 
at a constant density of 4 x 10 s per well. Every 7 days, cells were 



harvested and the total number was counted. The cells were 
then subcultured at the density described above. This cycle was 
repeated for five passages. Average PDs were calculated for each 
line and control at each passage according to the following 
formula: 

PD = (log 2) - l(logNt - logNO), where Nt = cell number at 
harvest and NO = cell number plated (4 x 10 5 ). The total PD level 
was calculated as the cumulative of average PDs for five passages. 
The PDT was calculated according to the following formula: 
PDT = I/PD, where I = interval (days) between the passages (7). 

Immunostaining 

Samples were processed as previously described (Ilic et al., 2007; 
Stephenson et al., 2012). The following antibodies were used: 
rabbit anti-K14, rabbit anti-filaggrin, rabbit anti-involucrin, and 
rabbit anti-loricrin (all from Covance); rabbit anti-laminin 
(DAKO); mouse anti-p63 (Santa Cruz Biotechnology); rabbit anti- 
fibronectin, mouse anti-collagen IV, mouse anti-collagen VII, and 
rabbit anti-desmocollin 1 (all from Sigma-Aldrich); and rabbit 
anti-p63, mouse anti-KlO, and rabbit anti-collagen I (all from 
Abeam). Secondary antibodies were obtained from Jackson Immu- 
noResearch and Life Technologies. 

Flow Cytometry with Confocal Imaging 

Cells were incubated with PE-conjugated rat anti-CD104 (integrin 
(34) antibody (BD PharMingen) for 1 hr, washed, fixed with 4% 
(v/v) paraformaldehyde, and permeabilized with 0.1% Triton 
X-100 in PBS. The samples were then incubated with fluorescein 
isothiocyanate (FITC)-conjugated mouse anti-K14 antibody 
(Abeam) for 1 hr, washed and incubated with Hoescht 33342 
(Life Technologies) for 5 min. To set the gating, fluorescence 
minus one (FMO) and single-stained compensation controls 
were also prepared. For acquisition, 1 x 10 6 cells/per sample 
were resuspended in 100 jxl of 1% BSA/PBS. The data were acquired 
on an ImageStream x analyzer and analyzed with IDEAS software 
(Amnis). 

Quantitative RT-PCR 

Data for quantitative PCR (qPCR) analysis were collected from 
three independently performed differentiation rounds for each 
of the cell lines. Total RNA extraction was performed at TO, Tl, 
T2, and T3 time points using the RNeasy kit (QIAGEN). For 
cDNA preparation, the QuantiTect Reverse Transcription Kit (- 
QIAGEN) was used. The primers are listed in Figure S7. KAPA 
SYBR FAST Universal 2x qPCR Master Mix (Kapa Biosystems) 
was used to set up the qPCR and run on a real-time PCR 7900HT 
system (Applied Biosystems). All reactions were run in triplicates. 
The data were collected and analyzed using GUSB and p-actin 



Figure 6. The Distribution of Epidermal Markers in HEEs Mirrors that in Human Skin 

(A) Immunostaining for three markers of keratinocyte terminal differentiation: filaggrin, Loricrin, and invoLucrin. Each of these markers is 
expressed at the appropriate site, denoting normal epidermal differentiation. 

(B-D) Immunostaining for K14 (basal layer) and K10 (suprabasal layers) (B), desmocollin 1 (suprabasal layers) (C), and p63 (basaL layer in 
general) (D) demonstrates normaL epithelial stratification in all HEE cultures. 

(E and F) Basal layer keratinocytes express laminin (E) and form BMZ-like structure (F). BM, basement membrane; Cy, cytoplasm; SB, 
stratum basale; TM, transwell membrane. Yellow arrowheads point to laminin-positive staining of BM. 
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Figure 7. HEE with a Functional Permeability Barrier 

(A) Lipid bilayer formation assessed with TEM. Arrows, Lamellar bodies; asterisks, Lipid biLayers. Upper row: normaL Lipid secretion between 
the SG and SCis seen in aLl cultures. MiddLe row: LBs are seen in the SG of aLL cultures (arrows). LB morphoLogy was normaL in aLl cultures, 
although it appeared slightly smudged in KCL034. Bottom row: lipid was successfully processed into lipid bilayers (asterisks) in all cultures. 
Like LB morphology, Lipid bilayer morphoLogy was normal in all cultures, although it seems slightly disrupted in KCL034. 

(B) ER Ca 2+ sequestration at day 7 of HEE culture in cells transfected with the ER-targeted Ca 2+ sensor D1ER. Data are presented as the mean 
of the intensity ratio (I) between the yellow-channel (higher Ca 2+ ) and blue-channeL (Lower Ca 2+ ) images (I(yfp)/I(cfp)) + SEM. Higher 
ratios denote higher Ca 2+ stores in the ER. n = 10-14 cells from two biologicaL replicates in each group. Significance was calculated using a 
one-way ANOVA. Distributions with p < 0.05 were assumed to be statistically different based on unpaired t tests between the populations. 

(C) The epidermal Ca 2+ gradient captured on TEM as electron-dense Ca 2+ deposits at day 14. Arrows indicate the direction from the basal 
layer toward the SC. Ca 2+ precipitates (black deposits) denote the presence of Ca 2+ in the tissue. Precipitates are seen in the viable SG but 
are absent from the SC, denoting a functional barrier to the passage of water and ions. 

(D) Permeability barrier integrity assessed by lanthanum perfusion. Lanthanum is visualized as electron-dense deposits in the extracellular 
spaces of the viable SG (arrowheads), demonstrating that lanthanum and, by extension, water and other small ions can pass between 
keratinocytes in this stratum. In contrast, lanthanum spreads along the base of the SC, but cannot penetrate further into the SC because a 
functioning lipid barrier is blocking its movement upward. ALl cultures demonstrated a functional permeability barrier. 
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or GAPDH as the reference genes. The results were analyzed and 
quantified using the ddCT method. 

Whole-Genome Gene-Expression Array 

Total RNA from T3 cell populations and control primary NHKs 
were reverse transcribed into cDNA using Superscript III Reverse 
Transcriptase (Life Technologies). In vitro transcription of 
biotin-labeled antisense RNA (aRNA) was performed with the 
TargetAmp-Nano labeling kit (Epicenter), and purification was 
performed with the RNeasy Mini Kit (QIAGEN). The quality and 
quantity of the resultant biotin-labeled aRNAs were assessed using 
a Bioanalyzer (Agilent) and Qubit Fluorometer (Life Technologies). 
A total amount of 750 ng of biotin-labeled aRNA in a 5 ul volume 
was then used for the HumanHT-12 Expression BeadChip 
whole-genome, gene-expression direct hybridization assay system 
(Illumina) according to the manufacturer's instructions and run 
on the iScan system (Illumina). 

HEEs 

A suspension of 0.25 x 10 6 neonatal NHKs or hESC/iPSC-derived 
keratinocytes was seeded on a CELLstart CTS (Invitrogen)- 
coated 12-well 0.4 um PET Millicell Hanging Cell Culture Insert 
(Millipore) in CnT.BM.l medium with CnT-07 supplement 
(CELLnTEC). Three days after seeding, the medium was switched 
to CnT-BM.3 medium with CnT-02-3DF with supplement 
(CELLnTEC). On day 4, the HEEs were air exposed by feeding 
the insert from bottom of the well only. HEEs were grown in 
a humidified (xl00% relative humidity) or dry (30%-50% 
relative humidity) incubator at 37°C (R.S., A.C., D.C., M.H., 
L.C., S.D. Peckenpacker, K. Park, Y. Uchida, K.R. Feingold, 
P.M. Elias, D.I., and T.M.M., unpublished data). The TEER was 
measured with an epithelial voltohmmeter (EVOM; World Preci- 
sion Instruments). 

ER Ca 2+ Sequestration 

Cells were transfected with the ER-targeted cameleon-based Ca 2+ 
sensor D1ER (a gift from Dr. Amy Palmer; Palmer and Tsien, 
2006) using the lipid-based TransIT transfection reagent (Mirus) 
at a 1:3 DNA/TrawsiT ratio. Transfection complexes were assembled 
30 min prior to cell plating by adding 1 ul D1ER (at a 1.3 ng/pl 
concentration) and 3 transfection reagent per 100 ul serum- 
free OptiMem. Cells were resuspended in low-calcium 154 CF 
medium supplemented with 0.03 mM Ca 2+ (GIBCO) at a density 
of 0.5 x 10 6 cells per 1.5 ml. Transfection complexes were added 
to the suspension at a 100 ul:l ml complex/cell suspension ratio 
and gently mixed by pipetting. The cells were then plated into 
12-well inserts precoated with CELLstart at a 0.5 x 10 6 cells per 
insert density and left with the transfection complexes overnight. 
The medium was then switched to CnT.BM.l medium with 
CnT-07 supplement (CELLnTEC) for 3 days. Dual-channel fluo- 
rescence measurements were recorded using a Zeiss LSM Meta 
confocal system (Zeiss). 

Transmission Electron Microscopy 

Samples were processed for transmission electron microscopy 
(TEM) as previously described (Ilic et al., 2007; Tu et al., 2012). 
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